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We theoretically study the skyrmion-number dependence of spin-transfer torque acting on mag¬ 
netic bubbles. The skymrion number of magnetic bubbles can take any integer value depending on 
the magnetic profile on its circumference and the size of the bubble. We find that the transverse 
motion of a bubble with respect to the charge current is greatly suppressed as the absolute value of 
skyrmion number departs from unity, whereas the longitudinal motion is less sensitive. 


In recent years, attention has been focusing on topo¬ 
logically nontrivial magnetic textures such as magnetic 
vortices [1] and skyrmions[2]. They exhibit rich physics 
stemming from their characteristic structures, which 
can be advantageous for technological applications [3]. 
Another interesting example among such topological 
textures is magnetic bubbles [4]; spot-like closed do¬ 
mains observed in ferromagnetic films with out-of-plane 
anisotropy, where the magnetization inside the bubble 
is oriented in the opposite direction to the one outside. 
Magnetic bubbles have a potential to play important 
roles in magnetic memory devices j4lll3]. 

Vortices, skyrmions and bubbles are quantified 
by a common topological quantity A/s, the so- 
called skyrmion number, which is defined by A/s = 
(l/ 47 r) J dxdy {m • dxTn x dyin), where m is the classi¬ 
cal unit vector in the direction of the local magnetization, 
and the integral is taken over the film sample. Whereas 
a vortex and a skyrmion carry A/s = ±1/2 and ±1, re¬ 
spectively, for a bubble A/s can take any integer value 
depending on the magnetic profile on its circumference 
and the size of the bubble. Dynamical response of a 
bubble to driving forces depends highly on its skyrmion 
number [HI [13]; a tantalizing prospect is that magnetic 
bubbles with different skyrmion numbers can provide 
a variety of new functionalities in device applications, 
which may not be obtained by skyrmions and vortices. 

In this work, we theoretically study the A/s-dependence 
of current-driven bubble motion. Micromagnetic simula¬ 
tions reveal that the transverse velocity of a bubble with 
respect to the current is strongly suppressed as |A/s| de¬ 
parts from unity, while the longitudinal motion is less 
sensitive. A collective-coordinate model (CCM), where 
the steady motion of bubble is assumed, provides good 
approximate solutions when |A/s| = 0 , 1 , and |A/s| ^ 1 . 

Let us begin by introducing the topological quanti¬ 
ties based on which magnetic bubbles can be classified 
[Fig. 1 ]; the winding number S counts how many full 
turns the magnetization on the perimeter of the bubble 
rotates, and its sign is determined by the sense of ro¬ 
tation. The polarity Q is defined to take ±1 when the 
magnetization inside the bubble points up, and —1 when 
it points down. The skyrmion number A/s is given by 
A/s = QS. Below we numerically examine the depen¬ 
dence of current-driven dynamics of a bubble on Q and 
S. The results of the simulation will be analyzed based 


on the CCM, where the mathematical expressions for the 
topological quantities are given. 

We assume that the magnetization obeys the Landau- 
Lifshitz-Gilbert equation; 
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where 7 is the gyromagnetic ratio, a and /3 are dimension¬ 
less parameters, iTeff is the effective magnetic field due 
to external, exchange, demagnetizing and anisotropy en¬ 
ergies, and u = —gyBPj/^^eMs with g the g-factor, gs 
the Bohr magneton, Ms the saturation magnetization, P 
the spin polarization of the conduction electrons, e the el¬ 
ementary charge, and j the charge current density. Here 
the charge current is assumed to flow in the x-direction. 

Eq. 0 is solved by the Object-Oriented Micromag¬ 
netic Framework simulator [TS] . where we divide a square 
thin film of dimensions 900 x 900 x 8 nm^ into 2x2x8 
nm^ unit cells, with the material parameters chosen to 



FIG. 1: Schematic of magnetic bubbles, where the black 
arrows indicate the magnetization. In the bottom figure five 
topologically different bubbles are shown [(a) S' = +1 and 
Q = +1, (b) S = —1 and Q = +1, (c) S' = +1 and Q = —1, 
(d) S = 0 and Q = +1, and (e) S ~ -\-R/A and Q = +1]. See 
the main text for the definitions of the symbols. 
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FIG. 2: Skyrmion-number dependence of the longitudinal 
and transverse velocities of bubble in the presence of charge 
current u — 100 m/s. Eqs. ^ are plotted by the open squares, 
while the results obtained by micromagnetic simulations are 
indicated by the red circles. The solid lines connecting the 
symbols are guides to the eye. See the main text for the 
material parameters used and the combination of Q and S 
for each A/s- The two calculations are in good agreement 
when |A/s| = 0,1 and |A/s| ~ R/A. 



FIG. 3: (a) Numerically obtained profile of a bubble with 

S' = +1 and Q = +1 in equilibrium. iiqHz = —2.7 mT 
is applied and ~ 110 nm. (b) Snapshots of the bubble 
motion over 2 ns after the current u = 100 m/s is turned on. 
(c) Trajectories of (X, Y) over 2 ns for three different current 
densities and Q — ±1. (X, T) = (0,0) is the initial position 
at t = 0, i.e., the center of the thin film. 


be typical for Co/Ni; 7 = 1.76 x 10^^ Hz/T, Ms = 
6.8 X 10^ A/m, the uniaxial anisotropy constant K = 

4 X 10^ J/m^, the exchange stiffness A = 10“^^ J/m, 

a = 0.03 and P = 0.01. A magnetic bubble is pre¬ 
pared at the center of the film in equilibrium, and an 
in-plane current is applied. We estimate the center-of- 
mass (X, y) of a bubble and its radius R by circular fit¬ 
ting to the numerically obtained magnetic profile [ 8 ]: X = 
Ei(i - Y = ^.(1 - 

TO|y /i(m^)/Ej(l - rn\)f i{m\), and R = y;.(l - 

+ -YYf2{m\)/Y.^{l - m*)/2(m*), 

where i denotes the unit-cell index, the weighing func¬ 
tion fi{m\) = 1 when |m* | < 0.99 but otherwise zero, 
and similarly / 2 (m*) = 1 when |m* | < 0.1 but otherwise 
zero. The bubble velocity is estimated from the displace¬ 
ment of (X, Y) divided by the time it takes. 

In Fig. 2, the results of the simulation is summarized; 
the bubble velocity is plotted by the red circles as a func¬ 
tion of A/s, where V\\ (W) is the longitudinal (perpen¬ 
dicular) velocity with respect to the current. (For the 
combination of Q and S employed for each A/s, see the 
discussion below.) It is clearly seen that \V±\ is greatly 
suppressed as |A/s| departs from unity, while V\\ is less 
sensitive to A/s- Below we will have a close look at the 
bubble dynamics at each A/s, and the results for |A/s I = O 5 
1, and 15 will be analyzed by the CCM. 

Fig. 3 (a) shows the equilibrium profile of a bubble with 

5 = -\-l and Q = + 1 , where the magnetic field ij^qHz = 
—2.7 mT is applied. R is estimated as 110 nm. A 


bubble with S = Q = —1 and the same radius can be 
obtained exploiting the reversed field iiqHz = +2.7 mT, 
see Fig. 1 (c) for a schematic. Fig. 3 (b) are snapshots 
of the time evolution of the bubble over 2 ns after the 
current u = 100 m/s is turned on. During the motion, the 
bubble sustains the circular shape and the magnetization 
profile shown in Fig. 3 (a). In Fig. 3 (c), the trajectory of 
(X, Y) is tracked for three different values of u and Q = 
±1 over 2 ns. The bubbles move with nearly constant 
velocities after the initial transient regime, and the travel 
distance is proportional to The sign change of Q leads 
to the change in the direction of the transverse motion. 
The results for |A/s| = 1 shown in Fig. 2 correspond to 
S = +1 and Q = ±1. The bubbles with S = —1 will be 
discussed later, where a qualitatively different behaviour 
than the bubbles with S = -\-l can be observed. 

Next, let us increase |A/s| to as large as 15. Fig. 4 (a) 
shows the magnetization configuration at the perimeter 
of a bubble in equilibrium where S = +15, Q = +1, 
IIqHz = —9.2 mT and R 110 nm. Shown in Fig. 4 (b) 
are the snapshots of the time evolution of this bubble over 
2 ns in the presence of current u = 100 m/s. The Bloch 
lines present in the domain wall region are so packed that 
the dynamics of the magnetization along the perimeter 
is suppressed enough to sustain the initial state’s profile 
shown in Fig. 4 (a) during the motion. In Fig. 4 (c), the 
trajectory of (X, Y) is plotted for the four topologically 
different bubbles all with |A/s| = 15 and the same ra¬ 
dius under the current u = 100 m/s; the direction of the 
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FIG. 4: (a) Numerically obtained profile a bubble in equi¬ 

librium, where S = +15, Q = +1, ijlqHz — —9.2 mT, and 
~ 110 nm. (b) Snapshots of the motion of the bubble over 
2 ns in the presence of current u = 100 m/s. (c) Trajectories 
of (X, Y) of the four bubbles with (S, Q) = (±15, ±1) driven 
by u = 100 m/s. (d) Charge current- and Q-dependences of 
the (X, y)-trajectories with S = +15. 


transverse motion is determined by the sign of A/s = QS. 
The linear dependence of the bubble velocity on \u\ is in¬ 
dicated in Fig. 4 (d), where the (X, y)-trajectories with 
S' = +15 and Q = ±1 are plotted for three different cur¬ 
rent densities. The results for S = +15 and Q = ±1 are 
shown in Fig. 2. 

A bubble with Q = +1 and S = A/s = 0, i.e., when the 
magnetic structure is topologically trivial [Fig. 1 (d)], is 
investigated with magnetic fields iiqHz = —11 mT and 
IJLqHx = 100 mT. The in-plane field is applied to lock the 
magnetization direction around the circumference in the 
x-direction. As shown in Fig. 2, is clearly suppressed 
compared to the other cases, while V±_ reaches zero. The 
small V±_ is not due to the actual translational motion 
of the bubble but due to its systematic deformation into 
an asymmetric ellipse accompanied by the shift of the 
center-of-mass. 

In the three cases discussed above, the bubble shape 
and the magnetization distribution along the perime¬ 
ter are rather rigid, motivating us to try to understand 
the results by a simple analytical model. Here we as¬ 
sume a perfectly cylindrical bubble with distribution of 
m = (sin 0 cos 0, sin 0 sin 0, cos 0) given, as schematically 
shown in Fig. 1, byg] 


0{r, X, z) = ±2 tan ^ exp 

4>{r,X,z) = Sx + (f>o, 


( 2 ) 

(3) 


where (r, x, 2 ;) is the cylindrical coordinate measured 
from the bubble center, A(^ R) is the domain wall width 
parameter, 0o is a constant. The topological quantities 
Q and S are defined by 



where f ds is the contour integral taken counterclockwise 
around the circumference of the bubble. It is straightfor¬ 
ward to prove A/s = QS. There are in general many 
possible ways of distributing the azimuthal angle 0 along 
the perimeter, and the linear dependence of 0 on x as¬ 
sumed in Eq. (§ is well satisfied only when |A/s| = 1 
[Fig. 1 (a)-(c)]. A/s = 0 with m at the perimeter aligned 
in one direction [Fig. 1 (d)], or |A/s| ^ R/A where the 
Bloch lines are packed so closely that the distance be¬ 
tween the adjacent Bloch lines is comparable to the do¬ 
main wall width [Fig. 1 (e)]. When a bubble contains 
a small number of Bloch lines, the ^-distribution is no 
longer as simple as Eq. (§. 

Let us employ (X, E) as the collective coordinate of 
the bubble dynamics. Assumed here is the steady motion 
of the bubble, where the bubble stays rigidly cylindrical 
with constant radius R during its motion, and the 0- 
distribution does not change with respect to the comov¬ 
ing coordinates. We Integrate Eq. Q over the sample 
volume [m [15] to obtain 


f V\\\ ^ u f \ 

\V^J + (aF)2 V Gr{a - /3) ; ’ 


(6) 


where m x iTeff = 0 has been assumed, and 


G = d^ATs, (7) 



The equation of motion of the same form with Eq. <11 has 
been known for a skvrmion [16111 7j . Eor a bubble, i) owing 
to the condition ^ A, which is usually not the case 
for a skyrmion, the analytical expression of F is accessible 
as in the second equality of Eq. (§, and ii) |A/s| is not 
restricted to 1, leading to the strong A/s-dependence of 
that enriches the bubble dynamics as already 
seen. Eq. (H is compared to the numerical results in 
Eig. 2 by the open symbols. The parameters are chosen 
to be consistent with the simulation at each A/s- The two 
calculations agree well at |A/s| = 0, 1, and 15. 

Lastly, we touch upon a couple of cases where the CCM 
is not a good approximation. A bubble with S = —1 
[Eig. 1 (b)] inevitably produces magnetic charges on its 
perimeter and thus is energetically unfavourable. This 
fact leads to relatively large shape distortion of the bub¬ 
ble during its motion, losing the legitimacy of using 
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FIG. 5: Snapshots of the time evolution of the bubble with 
aS = +3 and Q = +1 over 2 ns in the presence of current 
u = 100 m/s and the magnetic field Hz = —3.1 mT. 


two pairs of Bloch lines move along the circumference 
in the presence of charge current. The numerical results 
with S = +3, Q = ±1, and the Bloch-line distribution 
shown in Fig. 5 are compared with Eq. © in Fig. 2; the 
CCM completely fails to predict V±. We also observed 
that the bubble dynamics depends highly on the initial 
locations of the Bloch lines (not shown). The signal of 
the restoration of agreement between the CCM and the 
simulation is seen when |A/s| is increased to 7. We leave 
more systematic and complete investigations to future 
work. 

In conclusion, we presented analytical and numeri¬ 
cal studies on the current-driven bubble motion. We 
found that the transverse motion of the bubble with re¬ 
spect to the current is greatly suppressed as the bub¬ 
ble’s skyrmion number departs from unity. Our find¬ 
ings suggest the possibility to manipulate the dynamics 
of bubbles by their skyrmion number, which would lead 
to implementation of magnetic bubbles in wider range of 
applications. 


the CCM. By the micromagnetic simulation (not shown) 
with S = —1, Q = +1, fioHz = —2.2 mT, R 110 
nm and u = 100 m/s, we observed ^ 98.4 m/s and 
V± ^ —2.3 m/s; whereas agrees well with Eq. (1^, 
|Vj_| is about of an order smaller than the prediction oy 
the CCM. Shown in Eig. 5 is a case with As = +3. The 
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